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Abstract. This paper is concerned with the rigorous analysis of the zero electron mass 
limit of the full Navier-Stokes-Poisson. This system has been introduced in the literature 
by Anile and Pennisi (see [5]) in order to describe a hydrodynamic model for charge-carrier 
transport in semiconductor devices. The purpose of this paper is to prove rigorously zero 
electron mass limit in the framework of general ill prepared initial data. In this situation 
the velocity field and the electronic fields develop fast oscillations in time. The main idea 
we will use in this paper is a combination of formal asymptotic expansion and rigorous 
uniform estimates on the error terms. Finally we prove the strong convergence of the full 
Navier Stokes Poisson system towards the incompressible Navier Stokes equations. 



1. Introduction 

1.1. Model. In this paper we are concerned with the rigorous asymptotic analysis of the 
following scaled full Navier-Stokes-Poisson type system 

n t + div(nv) = 0, (1.1) 

A Ti TIV 

n[v t + v ■ Vu] + Vp me = div( — S(v)) + — qVV , (1.2) 

m e m e t p 

5 12 

d t (nk B T) + v ■ V(nk B T) + -(nk B T)divv + -divq me + -(\§(v) : Vv) 

o o o 

2 .1 1 ,, ,9 m e n , /m , 
= -nm e [- - —]\v\ 2 + -^k B (T - T), (1.3) 

AV = n-n. (1.4) 

Our goal is to investigate rigorously the incompressible type limits arising when m e tends to 
zero. This kind of incompressible limits are subject to various physical interpretations, the 
most immediate one, in some analogy with perturbative methods in quantum field theory, 
is related to the so called zero-mass limit in plasma physics (e.g. [16], [20], [4], [3], [9] ) or 
otherwise we may consider the quasineutral type limit arising when the Debye length is of 
the same order of the Mach Number [TTj. Other similar limits have been investigated by 

p], m- 

The system lll.l|) - (|1.4j) has been introduced in the literature by Anile and Pennisi (see [5]) 

in order to describe a hydrodynamic model for charge-carrier transport in semiconductor 

l 
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devices. Here is the list of notations within (jl.ip - (jl.4j) . 



A 



nksTT c 



nksT 



Qm e = -5 — k B nTr q VT + -(nk B T)[ }r q v, 

2 

SO) = Vv + (Vv) T - -divvl, 



m e 



where m e is the effective electron mass, ks the Boltzmann constant and t p , t w , r q are 
respectively the momentum, energy, total energy flow vector relaxation time. The aim 
of these models is to incorporate higher-order effects than those included in the standard 
drift-diffusion equations, in order to be able to describe high-field transport phenomena in 
semiconductors. These models includes the fundamental balance laws of particle density, 
momentum and energy for the charge carriers and are derived from the moment equations 
of the Boltzmann transport equation (BTE). 

In particular the right-hand sides of (|l.l)) - (|1.4j) representing the production of particles, 
momentum and energy, due to various interaction mechanisms (carrier-phonon, carrier- 
carrier and carrier-impurity collisions) are modeled by using relaxation type nonlinearities. 
Moreover in the Anile Pennisi (see [5]) derivation they ignore possible anisotropy of the 
stress tensor and they assume a Fourier-like constitutive law for the heat flux. The equations 
of momenta are closed by means of the principles of extended thermodynamics theory due 
to Mueller and Ruggeri [27J, the hydrodynamic systems then follows by the Maxwellian 
iteration method. This choice was motivated by them because of the success of the extended 
thermodynamics in obtaining the usual system provided by the Grad method for dilute 
gases. Such an approach allows to generate constitutive laws for the higher momenta of the 
distribution function, depending on the order of truncation, somehow independently from 
detailed microscopic assumptions imposed on the form of the distribution function. 

The mathematical analysis carried out in this paper gets its inspiration from the papers of 
[28] , [30J , [26] , [25] , [1] , where closed rigorous higher order expansions are carried out in order 
to manage the different oscillating components contributing to the formal asymptotics. As 
far it concerns the Low Mach number limit analysis there is a very large amount of relevant 
literature that we cannot report here in a reasonable way. Beyond the classical papers of 
[22] , [29] . [T7] . [18] . [3T] we mention the recent papers of [ij where the analysis presents 
certain similarity with our methods and the review papers of [2] . General references on the 
Navier Stokes equations are for instance the books of [231 El], [13] . 

According to the relation among all relaxation constants in Anile and Pennisi's derivation, 
see [5], 



5 T(jr 5 Tff > 



2t w > Tj 



1 



T, 



1 




we can take r CT = t p = r q = r, 



= 1. 
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Other physical constants, except the effective mass of electron m e , in system (ll.ip -( [L~4j) 
will not play any role in our analysis here, we can simply choose &b = Q = n = 1, and 
arrive at the following system 

n t + div(nw) = 0, 

11 n 

n[v t + v ■ Vu] H V(nT) = — div(nT§(») H W - nv, 



d t (nT) + v ■ V(nT) + -(nT)divu - --^-div(nTVT) + -(nTS(v) : Vu, 
3 3 2m e 3 

= — n\v\ + — n(T -T), 
AV = n-l. 

Moreover, for the convenience of analysis, let m e = e 2 , our system is changed into 

rit + v ■ Vn + nV • v = 0, 

Til 1 
u t + v • Vu + -^-Vn + -VT = — div(nT§(») + —VV - v, 

■rZryi £■£ £ ft *~ 



T t + -TV ■ v + v ■ VT 
3 

AV = n-l. 



JLdiv(nTVT) - l(TS(v) : V«) + |>| 2 + y(T - T), 



This hyperbolic system can be rewritten by using the following notations. Let U 
(n,v,T) T and F be the terms on the right hand side. 



U t + Y,Ajd Xj U = F, 



where the matrix Aj is defined by 



/ o 



A; 



nej 


-Tp ■ 







+ Vjl. 



\ §Te, 
It is easy to find that the multiplier to symmetrize the system is 





( 


T 








A = 







1 







\ 








3 

2e' 2 T 



such that AqAj 



/ t T 
T 



V 







\ 



1.2. Scaled system. A natural scaling of temperature T should be the one of the same 
order as electron mass. More precisely, T is of the order 0(e 2 ) such that Aq is uniformly 
positive. Now we use the new variable T = T/e 2 , and the same scaling for Tq = Tq/e 2 , 
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consequently we scale the electronic potential as V = V/e. Then the new system is 

n% + v ■ Vn + nV • v = 0, 

T 1 - 1 - 

v t + v-Vv + — Vn + VT = -div(nTS) + - W - v, 

n n e 

f t + -TV • v + v ■ VT = — div(nTVT) - -(TS : Vu) + -\v\ 2 + — (T - f ), 
3 6n v ; 3 V 7 3 1 1 2 v u ; ' 

eAV = n- 1. 

In the remaining discussion, for the convenience of our analysis, we will use the following 

Tl — 1 

group of notations, a £ = , u £ = v, T £ = T, Tq = To, ip £ = V, and the system with 

£ 

new notations reads 

of + ^divu £ + div(a 6 u £ ) = 

uf - -Vip £ + u £ ■ Vu £ = £ -!—-V(j £ - VT £ 

4 £ Y l + ea £ 

+ — !— -div((l + ea £ )T £ S(u £ )) - u £ , 
1 + ea £ 

Tl + ^T £ V -u £ + u £ - VT £ = — ^— - div((l + ea £ )T £ VT £ ) (S £ ) 

- ^T £ §(u £ ) : Vu £ ) + i|n £ | 2 + ^(T - T £ ) 

= a £ . 

2 

where E>(u £ ) = Vu £ + (Vu £ ) T - -divu e I. 

3 

The problem will be studied in a d dimensional periodic domain T d , with initial data 

a E \ t=0 = af(x), u £ \t=o = uf{x), T £ \ t=0 = If (x). (1.5) 

The main goal of this paper is to study the limit as e — > in fl^D - By a simple formal 
analysis, if we put e = in (Se_ ), we will get that n = 1, divu = and our limiting system 
will be an incompressible Navier-Stokes system with temperature of the form 

divv = 0, 

v t + v ■ Vv + Vn + VT - div(TS(u)) + v = 0, (1.6) 

T t + u • VT = ^div(TVT) - |(TS(«) : Vu) + ^| 2 . 

d 3 3 

We will give some detailed result on this system in Section 12.11 

The purpose of this paper is to prove rigorously this formal limit in the framework of 

general ill prepared initial data for the system fl^p . In this situation the velocity field and 

the electronic fields as e — ¥ develop fast oscillations in time. In order to describe them 

one has to introduce a linear oscillating limiting system that depends on the solution of 
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(11. 6p . We will give a detailed result on the solution of this fast oscillating system in Section 
1231 

The main idea we will use in this paper is a combination of formal asymptotic expansion 
and rigorous uniform estimates on the error terms. The ansatz of our solution is 

a £ (x,t) = (Ti f (x,-) + £[a 2 f(x,-)+<JE(x,t)], 

u £ (x,t) = v(x,t) + uif(x, -) + e[u 2 {(x, -) + UE(x,t)] , (1.7) 

T £ (x,t) = T(x,t)+e[T 2{ (x,^) + T E (x,t)]. 

Here we used notations for subscript "If" to represent the first order fast oscillation of 
the solution, "2f" to be the second order fast oscillation of the solution, and "E" to be the 
error terms. 

In the above expansion, (v,T) is the solution of limiting incompressible system (|1.6[) . 
(o"if,uif) is the solution of first order oscillation system and (o"if , u 2 { , T 2 f) is the solution of 
second order oscillation system. The exact form and the solvability of these systems will 
be listed in the coming discussions. 

Since we have the relation a £ = Atp £ , we have similar representation for electronic 
potential 

V^ £ (x,t) = V^2f(x,^) + e[V^(x,-) + V^E(x,t)]. 
For consistence, we also write down initial data in the same form. 

(Tf(x) = Alpl{x)+£uf{x), 

uf(x) = vi(x) + Qui(x) + euf \x) , 
T{{x) = T^+eTfix), 

where Q is the Leray's projector on the space gradient of vector field u G L 2 (T d ) and is 
defined as follows, 

Qu = VA _1 divu, Pu = (I - Q)u, divPn = 0. 

So, our task is first to find the limiting system, which will be the incompressible system 
(II. 6p with solution (v,T, II) and a fast oscillation system with solution (uif,<7if = A^if). 
Under suitable assumptions on initial data v\(x), Tj(x), divv\ = 0, we can prove that fjl .6|) 
is solvable up to time r*. Moreover, for any s > f + 1, r G (0,r*), (v, T) G L°°([0, t];H s+3 ) 
and n G L°°([0, r]; H s+3 ), d t U G L°°([0,t]; H s+2 ). We can also prove that the solution of 
fast oscillation system satisfies (a 1{ ,u 1{ ) G L°°([0,t]; H s+3 ), V^if G L°°([0, r]; H s+3 ). 
Now we are ready to state the main result of this paper. 
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1.3. Main result. 

Theorem 1.1. Let be s > ^ + 1. Assume that the initial conditions (jl .51) satisfy 

||A^|| H . +1 + ||Qui|| H .+i + \\eaf\\ H s + \\uf \\ H . + \\Tf\\ H s < C. 

Then, there exists an eq > 0, such that for all e < eq, problem (\S £ \j fjl .5f) has a unique 
classical solution (a £ ,u £ ,T £ ) G L°°([0, r]; H"), u £ ,T £ G L°°([0, r]; H s ) n L 2 ([0, r]; 
and 

sup - trif)||H. + ||u £ - u - u xt \\ H . + \\T £ - T\\ H . + ||V^ e - V^ f || H .] < Ce. 

0<t<T 

\\U £ -V- Ulf ||i2([ 0)T ];Ha+l) + ||T e - r|| i2 ([ 0iT ].^ a +i) < Ce. 

Our paper is arranged as follows. In section 2, we will discuss the limiting incompressible 
system and the leading order oscillating system and give the results on their existence. In 
section 3, we use asymptotic expansion to find a second order oscillating system and give 
the derivation of the error system. In section 4, a detailed proof on existence of solution 
to the error system by using energy method will be given. In the appendix, we give a 
motivation on how to get our leading order oscillation system. 

2. Limiting system 

2.1. Limiting incompressible system. In this section we give some more details on our 
limiting incompressible Navier-Stokes system, namely 

divv = 0, 

v t + v ■ \7v + VII + VT - div(TS(u)) + v = 0, (2.1) 
T t + v • VT = ^div(TVT) - -(TS(«) : Vu) + ^M 2 . 

with initial data 

v\ t=0 = v l (x), dwv l = 0, T\ t=Q = T I (x)>T h >0. (2.2) 
For the system (|2.ip the following existence results for smooth solutions holds. 

Theorem 2.1. Ifv h Ti G H s+3 and Ti > T L > 0, then there exists r* > such that (|2~1) 
(|2.2p has a unique solution such that for any r < r* the following holds 

SUp [|K-,*)||tf»+3 + ||T(-,t)|| HS +3] + \\v,T\\ L 2n ]H s+4) 

0<t<T 

<C(t,T Jj )(\\v 1 \\ Hs+ s + (2.3) 

SUp [\\U(;t)\\ HS+3 + \\dtU(;t)\\ H s +2 ]<C(T,T L ). (2.4) 

0<i<r 
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The proof can be obtained by extending the result by Kato [21] for the incompressible 
Navier Stokes system to the full system and by using the parabolic theory for the tempera- 
ture equation. For completeness we mention here the following existence results regarding 
the coupling of the incompressible Navier Stokes equations and the temperature balance 
equation, [ID]. [8]. [6], [7]. 

2.2. Leading order oscillation system. By following the same line of arguments as 
Masmoudi's |25] and Schochet's [28] we describe the main fast oscillating system as follows, 

2d t Vq = Q{-(Vg-V)u-(u-V)Vg + vAg + div[TS(Vg)]}-Vg, (2.5) 
2d t V(f) = Q{-(V<p-V)v-(v- V)V0 + vA0 + div[T§(V</>)]} - V0. (2.6) 

with initial data 

Vq\t=o = Qui(x), V<P\ t=0 = V^i(x), (2.7) 

where v is a divergence free vector field. The details on the construction of such a system 
can be found in the appendix. 

Then, our leading order fast oscillating system can be obtained from (|2.5|) and (|2.6p by 
setting 

where for any v, e € L 2 (T d ; M. d ) with divv = dive = 0, L is defined as follows 

and L ( M = (° -' )( V " V 
Accordingly the system for the fast oscillating vector fields (o"if,wif) is given by 

cVif + -div{(V^ lf • V)u + {v ■ V)V^if-w lf - div(TS(VV>if))} (2.8) 

+ o"if + -divnif = 0, 

e 

dtuu + -Q{(u 1{ ■ V)v + {v • V)ni f - vV ■ u 1{ - div(TS(u lf ))} + u u - -V^u = 0, (2.9) 
and the equation for S/ip ai{ is 

$Wif + ^Q{(V^ lf • V)« + (u • V)V^if - uA^i f - div(TS(V^i f ))} + V^if + -«if = 0. 
The initial data of this oscillation system are 

VVrf|t=o = V^i, uif|t=o = Qui, crif = A^i. (2.10) 
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We get the existence of the leading order oscillation system immediately from Lemma 
16. II in the appendix. 

Lemma 2.1. The Cauchy problem (12.9p(12.8p ( 12.10p has a unique solution (<7if , u%{ , V^if) 
which satisfies the following uniform in e estimates 

SU P [lklfllif+3 + |klf||#s+3 + HV^lfHiJ-s+4] < C(t). (2.11) 

0<t<T 



3. Expansion of the system (S< 



3.1. Second order fast oscillation system. By direct calculations from the expan- 
sion (jl.7p . and by comparing the system with a combination of its incompressible sys- 
tem and the leading order oscillation, we can get that the second order fast oscillations 
(cr 2 f , U2f , T 2 f , VV^f ) satisfy the following system, 

d s a 2 { + divu 2f = F 2f ', 

a s n 2 f-VV 2 f = Fl\ 

d s T 2f = F| f , (3.1) 

A^2f = 0"2f, 



with initial data 

C2f|s=0 = U2{\s=0 = 72f| s=0 =0, (3.2) 

where the right hands of the system is given by 

Ff = idiv{(V^if •V)v + (wV)Wif-wif-div(TS(V^if))} 

+cr lf - div(crif(v + uif)), 
Fu = \q{{uu ■ V)v + (v ■ V)«if - vV ■ u 1{ - div(TS(« af ))} 

— u\fV(v + nif) — v ■ Vtiif + div(TS(«if)), 

2 2 
F| f = --TVuif - uifVT - -TS(u)V«if 

o o 

-|rS(uif) : V(« + u u ) - ^M 2 + + uu\ 2 . 

For the linear system (|3.ip . one easily proves the following result 

Lemma 3.1. Problem (|3.ip(j3.2p has a unique solution which satisfies the following uniform 
in e estimates 

sup [|k 2 f(-,s)||^+i + \\u2t(;s)\\ H s+i + \\Vip 2 {(;s)\\ H s+2 + ||T 2f (•, s)|| H »+i] < C(r). (3.3) 

0<s<oo 
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3.2. 0(e) Approximation system. Summing up the limiting incompressible system (12. ip . 
the first order oscillation system f)2.8j) ()2.9|) and the second order oscillation system (|3. 1 1) . 
we have that the 0{e) correction of our scaled system ( S £ ) is 

<9t(o"if + ea 2 {) + -div(-uif + eu 2 {) + div(<7if (v + uu)) = 0, 

d t (v + nif + eu 2 f) - -V^if + £ip2i) + (v + nif) 

+ (v + u 1{ ) ■ V(« + u 1{ ) - div(TS(t> + u lf )) + Vn + VT = 0, (S lf ) 

d t (T + eT 2f ) + \TV ■ (v + uif) + (v + u lf ) • VT 

+ ^TS(v + uif) : V(« + uif) - ^div(TVT) - hv + u lf \ 2 = 0. 

3 O 3 

We will use the following notations for the approximation of the solution. 

<7app( X ' *) = a lf + ecr 2f f ^ 

-"■a PP (M) = v(x,t) +uif (x, ^\ +£U 2 t ( ' ■ - ) (3.1) 

T app (x,t) = T(x, i) + eT 2{ fx, ~ 
The initial data of our approximation problem is 

o" app |t=o = A^i(x), u app | i=0 = vi(x) + Qu\(x), T app |t =0 = Ti{x). (3.5) 

Moreover, by the results of Theorem 12.11 Lemma 12.11 and Lemma 13. 1\ we have the 
following uniform estimates for this approximation solution 

SUP [||<7app||H»+i + ||«app||tfs+l + || V0app II-EP+ 2 + II 

T app ||//s+i < C(r). (3.6) 

0<t<T 

T ap p > T L /2 > 0. 

3.3. System for error terms. Now, by comparing O(e) correction system ( S\ f ) with our 
scaled system |52), we get the following system for the error terms (cte, u e , Te, ~4>e), 

l + ecr £ ^ E , r E 



d t (a E + An) + uapp • V(<r E + An) + — - — div^E = G„ + 

_£7^ v((te + An) + vr E + u E - - 

1 + ea e e 

iE 



d t u E + n app • Vu E + 1 , a !L v (^E + An) + VT E + u E - - V(^e + n) (S L 



-div(T w S(u E )) = G£ + F*, 
d t T E + ^T app V • u E + u app • VT E - ^div(T app VT E ) = Gf + if , 
with initial data 

^ E \t=o = &\ , u E \ t =o = uf, T E \ t=0 = Tf. (3.7) 
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The right hand sides of the error system (\Se\ ) are given by 
= -V(<7 app + ea E ) ■ u E - cr E divu ap p, 
Gf = -u E • V(n app + eu E ) + div(T E S(u app + eu E )) + gTa PP v(a E + An) 

1 + £(T t 

T £ V<7 £ 

:V/ + - -T £ S(n e ), 



1 + eo e 1 + ea £ ' 

G# = ~\t e V ■{u &w + eu E )-u E - V(T app + er E ) + ^div(r E V(T app + er E )) 



2 2 
■-e 2 T E S(n E ) : Vu E - -eT E S(it E ) : Vu app 

2 2 

--eT E S(u app ) : Vu E - -eT app S(u E ) : Vu E 

2 2 2 

-T app S(u E ) : u app - -T app S(n app ) : Vu E - -T E S(-u app ) : Vu app 

e, ,o e 2 „ 5Vo- e 2^ 

"'app 



+^e| 2 - ~T E + „„ — ^ T e VT £ + ^, ip! , • » K . 

and 



2 6(l + ecr £ ) 3 



F CT = -div(cr 2 fu a p P ) - div(«7if« 2 f) + d t ALT + -u app VAn. 

F„ = -u 2 f - VT 2f - u 2 f • Vu app - (v + ui/) • Vn 2f + div(T 2 fS(u app )) + div(TS(n 2f )). 
2 2 

= --^V • u 2f - -T 2f V • n app - u 2f • VT app - + uh) • VT 2f 

2 2 2 

--r 2f S(ti app ) : Vu app - -TS(n 2f ) : Vii app - -TS(v + u lf ) : Vu 2f 

5 5 2 e 

+-div(T app VT 2 f) + -div(T 2f VT) + -ti app • u 2f + (v + iiif)it 2 f + ^ ( T o - r api 

Now, the last step is to prove the following existence result for the error terms. 



Proposition 3.1. Let s > f + 1. Ifaf,uf,Tf G H s , then flffjl) (|3.7j) /ias a unique solution 
such that 

sup ||ea E ,u E ,T E , V^ e ||h s + II^e^bII^qo i Hs +i) < (7. (3.8) 

0<t<r 

The proof of this proposition will be done in the next section. 

4. Existence of solution to the system <\Se\) of the error terms 



In this section, we focus on solving the system of error terms by using energy method. 
The following Moser's type inequality will be used several times in the coming estimates. 
For completeness, we list it here as a proposition. 

Proposition 4.1. The following facts hold 
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(1) Let s > 0, f, g G H s (T d ) n L°°(T d ), and a a multi-index with \a\ < s. Then, for 
some constant c s > 0, 

\\D a (fg)\\ < cdl/Hooll^llo + |b||oo||-D7l|o). (4.1) 

(2) Let s > 1, f G H s (T d ) with Df G L°°(T d ), g G H s ~ l (J d ) n L°°(T d ), and \a\ < s. 
Then, for some constant c s > 0, 



ll^jbllo^c.ciP/iuii^-^iio + lblUII^VIIo)- 

where [D a , f]g = D a (fg) - fD a g. 



(4.2) 



We rewrite the system QS^D into a symmetrizable hyperbolic system. We introduce some 
notations to shorten our formula. 





( 


o-e + An \ 




( 





e e 3 





U E = 






, Af = 


i- 







e j 




\ 


T E J 




\ 











Then we can choose the symmetrizer of the system to be 

\ 

and Af = A E A~ 



A 



\ 



T 











/ 











3 

2Xapp 



V 



£ ^~app , 

l+£cr^ C J 

o 



en. 



app/J 



J 2T a , 



;(^app)j / 



Therefore, the error system ( Se ) is written down into a symmetric hyperbolic system 
with viscosity 



d t U E + Afd Xj U E - VUe = —-J + H 



where 



VUe 







(4.3) 



\ 



-u E + div(r app S(u, E )) 
V §div(T app VT E ) J 



J 



\ 

-V(Ve + n) 
/ 



and 



H E = H§ + H§ 











V 


F U E 









The symmetric version of our system (Se) is 



A%d t U E + Afd Xj U E - A$VUe = —A$J + A E H E , 



(4.4) 



In the following we will use the operator 



OJ = A E 8 t U + Afd Xj U - A E VU, 
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and, for any multi- index a, we will use the notations 

a a = D a a E , ip a = IT E + II), u a = D a u E , T a = D a T E and U a 



with the following relation a a = Aip a . 

Now for any multi-index a, < \a\ < s, s > i + 1, by applying the differential operator 
D a on the hyperbolic system (I4.3P and by symmetrizing it, we have 



CU a = A$d t U a + Afd x .U a - A^VU a = --A$D a J + TZ a (4.5) 



TZ a = A$D a H E + A^[Af,D a ]d X] U E -A^[V,D a ]U E , (4.6) 

where we used the commutator notation [•, •]. 

By applying the energy method for the system ()4.5f) . we obtain 

(CU a , U a ) = ~\ j • Ua + (K a , U a ). (4.7) 

where (f,g) = J f 9 and also 
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J (I^qI 2 + \u a \ 2 + \T a \ 2 ), and ||f7 E ||a,e = ^2 W Ua 



\a\=0 



Now we divide the estimates into three steps. In the remaining estimates, we will denote 
C to be constants depending only on the following quantities, 

SUp ||f/ E || Sl e, Iklfj 02f lli°°([0,r];H a + 1 )! II v i U U > u 2f \\l°°{[0,t];H s + 1 ) > 

0<t<T 

||r,T 2 f |U°o([0,r];H s + 1 )) l|n||Loo([o,r];// 3 + 3 )> ||^fcH|| L oo ([Q )T ] ; /pr+2) , Tl . 

Step. 1. The left hand side for linear operator C in (|4.7|) yields to 

Id f , £ Z~app 2,1 1 2 • ^ 



(cu a ,Ua) = -- /( ,, , dp :, 9 <+K»r+ 



app 

' 2 

'- a I 



+ y |-"a| 2 + J T app §(ii a ) : Vm q + - J |VT Q 

+ I / ^— TaVT -PP VTa ~ / ( 9 ^0 U a ) Ua-^J (d Xj AfU a ) ■ I7 a , 



-app 

where the last three terms can be estimated by 

eWJIe + 7 / |vtj 2 . 
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Hence we get 



+ j K| 2 + y J \Vu a \ 2 + J \VT a \ 2 - c\\u a \\l (4i 



Step. 2. The singular term on the right hand side of (|4.7p can be handled by using the 
first equation of (|4.5|) . i.e. 



8 t A^ a + D a (u app • VA^o) + ^dwu a + a £ dWu a = [1+^1, D a ]divu E + D a (G^ + F?). 



Now the singular term in the energy estimate f|4.7j) is 
If If 

/ Vlp a ■ U a = - / ^ a dlVU a 

£ J £ J 

= - J $ a dtM>a ~ j ^aD a (u app ■ VAV'o) - J ip a a £ divu a 

+ J M^-^,D a ]dw U v + J i> a D a {GZ + FZ) (4.9) 



1 d 
2~dt 



j \Vip a \ 2 - J V'aWapp • VA^ a - J ip a [D a , u app ] VA^ - J ip a cr £ divu c 



where the fist term on the right hand side contributes in the energy. 

The last four terms in (|4.9p are easily estimated by the energy as we will do later, while 
the second term is different. Since we have only ea a = eAip a in the energy, there is no hope 
to control it just directly by Aip a and Vip a . Our idea here is to move one of the derivatives 
onto M a pp which is known and W 2 '°° controllable. The reason we can succeed in doing this 
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is that this term is quadratic in tp a . More precisely, 

- J V^app • S/Aip a = J S/ip a ■ U app A7p a + J ^adivn ap p Aip a 

= J diip a ul pp djjip a + J V'Qdivu a ppAV' Q 

= ^2 ®i3i>otUlppdjll>a + ^2 / dilpadjlllppdjlpa - / V(^ a divn app ) Wo 

i,j=l J i,j=l J J 

= \ it j <PP d i( d ^) 2 + / (W«Vu a pp) • VV> a 



VV' a divu a ppVV' Q : - / ip a Vdivu app Vip a 

= ~\ j divUappjVV'al 2 + J (VV'aVUapp) • Vlp a 

- J VV> Q diviia PP VV> Q - J ^ a Vdivn app V^ a . 
Thus this term can be estimated by 

- / V^app • VAl/) a < ||ttapp|| 

The estimate for all the terms, but the last one, on the right hand side of (I4.9P can be 
done by integral by parts and Moser's type inequalities Proposition 14.11 i.e. 



-J 1p a [D a ,U aiPp ]VAlpo < ||« a pp|lw 2 .°°IIV'o|lH3, 

J jp a a £ 6xvu a < C(||e«7 ||H* + ll^allia) + ^ J |Vu 



and 



lpa[ ,D JdlVU E 



< C(ll0a|ll,2 + Iko-ollff* + \\ue\\hs). 

The last term in (|4.9|) can be estimated by using integral by parts, Holder's inequality 
and Moser's type inequalities, Proposition 14. H 

< C(\\V^ \\ 2 H s + \\ea f H s + \\u E \\ 2 HS ). 
Then our estimates on the singular term in ()4.7p is 

Vip a -u a <~ J divnapp|VV Q | 2 + C(||C/ E ||f iC + ||VVo||^) + ^ J |Vn Q | 2 . (4.10) 
Step. 3. Finally, we estimate the remainder terms (1Z a , U a ) in (|4.7p . 
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The commutator terms are controlled by using Moser type inequality, Proposition l4.lt 



J A^[Af,D a ]d X] U E U a -j A*[V,D a ]U E U a 



<C(||ue||^+i + ||Te||h s+ i)||[/ q 



The last error term is 



J A§D a H E -U a = J A%D a H% -U a + J A§D a H^ ■ U a . 

Since we have the viscosity and the heat diffusion terms in our system, we can move one 
derivative from the nonlinear term to u a and T a . While for density o a , we couldn't do this. 

Now we know that all Hp is given by our approximation function which are known. We 
can simply estimate the second term above by 

J A^D a Rp ■ U a 

< C{\\D a F E \\ L 2\\ea \\ H s + \\D a ^ F E \\ L 2\\Vu a \\ L 2 + \\D a - x F^\\ L 2 \\VT \\ H s) 

< C + C\\eao\^ H s + ri\\V u a \^ L 2 + ^||VTo||^- s , for small 77. 

The estimates on Hq, G e and are relatively easier by using Moser 's type inequalities, 
Proposition 14. 1[ Just notice that whenever we have a D a+1 on ue or Te, there is an e in 
front, which give us the possibility to control them by using the viscosity and the heat 
diffusion terms. The result is 



J D a G*u a + f D a G\ 



2T f 



app 



< 



I 7-ia— 1/~<E II 1 



V« Q || L2 + \\B a G^\\ L 2\ 



2Z 



'allL 2 



app 



< C + C\\U E \\l e + \\eVu a \\ 2 L 2 + \\eVT a \\ 2 L2 . 

But we have to be more careful with the term D a G E , since we could not move the 
derivative to eo a since there is Vcte in which will give us D a+l o E and we will not close 
our estimates with this. Our method here is similar to the one we used to deal with singular 
term (|4.9p . More precisely, we have 



/ 



D a Gt 



e 2 T., 



app 



(1 + so") 2 



Or 



< 



< 



< 



D a {eVo E • « E ) 



e 2 T { 



app 



Or 



eVo a ■ u E 



e 2 Z 



app 



V(T Q -u e - 



(l + eo £ ) 2 
e 2 T f 



+ 



-D Q (V(j a pp • u E + <7 E divu 

e 2 Z 



e 2 T f 



a|>|) 



O r 



e[D a ,u E ]Vo E 



-app 



(l + eo £ ) 2 



+ C\\U E \\ 



2 

s.e 



app 



1 +£O e ) 2 
£ 2 T 



+ C\\U E 



|2 

ls,e 



-app 



-K| div^ E(1 + £(7£)2 



+ c\\u E \\l e <c\\u E \\l l 
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According to the previous estimates, we arrive at 

(K a ,U a ) <C + C\\U E \\l c + v\\Vu a \\l2+ri\\VT \\ 2 H s, for small r, > 0. (4.11) 
Finally, from (I4.7P , (|4.8p , (14.10P and (14. lip , we end up with the energy estimate 

d f , £ Ta rm O i , O 3 



+ J \u a \ 2 + J |Vn Q | 2 + J \VT a \ 2 <C(\\U E \\l e +\\V^ \\h) + C. 
Taking the summation for all < \a\ < s and by using Gronwall's inequality, we have 

sup (||Be||L + WVMh) + / (IMIfP+i + IITeII^+i) < Ce Cr . (4.12) 

0<t<T JO 

5. Proof of the Main Result Theorem 11.11 

By using the asymptotic expansion ()1.T|) . the Theorem 1 2. H and the Lemma [2.1l we get the 
existence and uniqueness of classical solutions to the initial value problem for the system 



([5^1), f|1.5j) and the solution satisfies 



sup [\\e(a £ - a u )\\ H s + \\u £ -v- Uli \\ H s + \\T £ - T\\ H . + ||V^ e - V^ lf ||^] < Ce. 

0<t<T 

\\u £ - v - u X i ||i2([ 0jr ] ;ff 3+i) + \\T e - T|| L 2([ 0iT ].^+i) < Ce. 

6. Appendix: Formal derivation of leading order oscillation 

We follow the ideas of Masmoudi's [25] and Schochet's [28] in order to deal with the 
leading order oscillation system. We start from ( [571 , for simplicity, we drop all e in the 
superscripts. 

First, it is reasonable to rewrite the mass conservation equation by using the electronic 
field 

dtVijj H — u = —{<?u) = —uAi/j. 

Since there are no oscillation for the emperature equation, the leading order oscillations 
have to come from the equations from mass and velocity, namely 

*( " )+-(° " )J — V " + F ("> + G ) ,6.1) 



where 



1 



F(u) = F(ea, T, u) = div (1 + ea)T(Vu + (Vu) 1 - -divul) 



2 



1 + ea 

G = G(sa,T) = — Vd-VT. 

1 + ea 



3 



+ u, 
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6.1. Mapping L and e Lt . Let L be a linear mapping from L 2 (T d ;M. 2d ) to itself, which is 
defined in the following way. 

For any ti,e£ L 2 (T d ; M. d ) with divu = dive = 0, we set 

£ (0 =£ (°) =0, (6 - 2> 

Since for any (u, E) T G L 2 (T a! ; M. 2d ), it is well known that the Hodge decomposition yields 

to 

where v,e are divergence free and q,4> G ff 1 (T ci ; K d ), L is well defined by the identities 

412]), 463}. 

The eigenvalues of L are 0, i and — i, the corresponding eigenspaces are given by 

E = |(o)'(°)' Where V,eG L2 ( jd ' and d[vv = dive = °| 
Et = | y j , where g G iJ 1 (T d ;]R d )| 

= | ^ j , where g G i7 1 (T a! ;M a! )| . 

For any (u,E) T G L 2 (T ci ;]R d ) with -u = w + Vg, E = e + V<^>, the projection operators 
according to L are given by 



Pn 



P 



£ J I e 



u \ 1 / Vg + iV^ 

£ J 2 I -iVg + 



P. 

Then (/ — Po)e rL is defined by 
(/-Po)e rL 



u \ II Vq-iV</> 

El 2 I iVg + 



u \ / cos rVg — sin rV0 \ / cost — sinr \ / Vg 
E 1 / \ cos rV0 + sin rVg / \ sin r cos r / I 
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6.2. Asymptotic expansion by fast and slow time. Assume that we have the expan- 
sion of solution, let r = | 

u(x,t) = u (0) (x,t,r) +eu (1) (x,t,r) + e 2 u {2 \x,t,T) H , 

Vif>(x,t) = Vi/j {0) (x,t,T) +eVV' (1) (x,t,r) + e 2 V^ (2) (x,t,r) + ••• , 

T(x,t) = T {0 \x,t) + eT w (x,t,T) + e 2 T i2 \x,t,T) + --- , 

cr(x,i) = Aip(x,t) = A^ \x,t,r) +eAV' (1) (x,t,r) + e 2 A^ (2) (x,i,r) H , 

By putting this ansatz into our system and comparing the orders of e, we have the following 
systems. 

The O(i) order terms give 

d rU (o) _ vv> {0) = 0, 

d r VV> {0) + W°> = 0. 

The O(l) order terms give 

dtU (o) +drU (i)_ Vljj (i) = - u W. V u^ + F(0,u^,T^) + G(0,T^) 
W (0) + 5 T VV {1) + = -u^A^°\ 
By using operator L we defined before, we are able to rewrite this system into 



u 



(i) 



,-£t 



(i) 



+e 



-Lr 



uW(x,t,0) 
VV> {1) (M,0) 



_a^(o) _ u (o) . Vn (o) + F(0) u( o) ) T (o) ) + G ( 0) r(o) ) 



gLs. 



To have the expansion meaningful, the first order term should be , Vip^ Crasr-^oo. 
Since it is obvious that the initial data term 



-Lt 



uW(i,t,0) 

vv- (1) M,o) 



o, 



then it is necessary to check that 



1 

T JO 



-d t uW-uW-VuM + F(0,u( 0) ,TW) + G(0,TM) \ 
_a t VV {0) -u(°)A^°> 



) 



ds^O, as t — y 0.(6.4) 



From the eigenvalue analysis of operator L, we know that for any given {u^°\Vip^), 
there exists (Vg, V0) such that 



( 



u (0) 



o 



and (J - P„) 



(0) 



I w (0) 
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In other words we have 

n (o) \ ^ / v \ _ tL / Vg \ _ / v + cosfVg + sin|V0 \ 
V^(°) ) ~ \ J** ' \ V<f> ) ~ \ cos f - sin f Vg J ' 

By applying operator Pq on both sides of (|6.4p . we will have the incompressible limit. 
Applying operator / — Pq on both sides of (|6.4p . we will get our leading order fast 
oscillation system, which will be shown in the following calculation. 
We will calculate the following limit 

lim - f dr(I - P K L ( Qh ) 

lim 1 [ S dr ( C ° S T ® Ix ~~ T< ^ 2 \ 
s^+oo s J y cos tQI 2 + sin tQI\ J 

J_ / cosrQ/i -sinrQ/ 2 \ 

27r 7 I cos tQI 2 + sin tQI\ J ' 

where 

Ji = — (v + cosr Vg + sin rV</>) • V(i> + cost Vg + sin rV^) 

+ cos rVg + sin rV</>) + G 
-^2 = — + cosr Vg + sin rV0) • (cost A</> — sin rAg). 

Thus the desired oscillation equations are 

28 t Vq = Q{-(Vq-V)v-(vV)Vq + vAq + div[T§(Vq)]}-Vq, (6.5) 
2d t V(/) = Q{-(V<f>-V)v-(v- V)V0 + i;A0 + div[TS(V</))]}- V<f>. (6.6) 

with initial data 

Vg| t=0 = Qui(x), V^|t= = Vifnix), (6.7) 

where ^i(x) satisfies A^i(x) = a\{x). 

It is straight forward to prove the following result by using energy estimates. 

Lemma 6.1. For any fixed t > 0, Vs > | + 1. Given v,T 6 L°°([0, r]; # s ), T > T L > 
and Qtii(a;), V^i(x) € -£P. Problem (|6.5p(|6.6p(|6.7p has a unique solution (Vg, V</>) suc/i 

sup ||Vg(-,t), V0(-,t)||^ + ||V?,V0|| r 2 ([O)T] . H . +1) < G(t)||Q«i, WiIIh- (6.8) 

0<t<T 
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